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Abstract 

The passage from angular to spatial correlations, in the case of spatial clustering 
length depending on the average distance between nearby objects is studied. We 
show that, in a number of cases, the scaling law of angular correlation amplitudes, 
which holds for constant spatial clustering length, is still true also for a luminosity 
dependent spatial correlation. If the Limber equation is then naively used to ob- 
tain 'the' spatial clustering length from the angular function amplitude, a quantity 
close to the average object separation is obtained. The case of cluster clustering is 
explicitly considered. 
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1 Introduction 

Initial studies on clustering properties over large scales were performed using 
2-dimensional samples (Totsuji & Kihara 1969, Peebles 1980; see also Sharp, 
Bonometto & Lucchin 1984). Spatial properties were then deduced from an- 
gular ones using the Limber equation (Limber 1953). However, even in the 
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present epoch of large redshift samples, the technique used to pass from angu- 
lar to spatial statistical properties is far from becoming obsolete. For instance, 
cluster samples are available and in preparation, selected with a fixed apparent 
magnitude limit. Even though the typical redshift of each cluster is given, the 
different intrinsical luminosity limit at different distances forces to cut off some 
information, to extract volume and intrinsical magnitude limited subsamples, 
where the spatial 2-point function can be directly estimated. Working out 
the angular function and using the Limber equation to extract spatial proper- 
ties, might then seem a more convenient alternative approach. This approach 
was followed, e.g., by Dalton et al (1997), when discussing the clustering of 
clusters in catalogues extracted from the APM Galaxy Survey. In the final 
section we shall comment on the results of their analysis. The point is that 
this approach meets a different obstacle, as the cluster correlation length Tc 
may depend on cluster luminosities L (see Bahcall & West 1988 for a review 
of Abell cluster clustering data; see Nichol et al 1992, Dalton et al 1994, Croft 
et al 1997, for APM cluster clustering data; the whole set of data has been 
recently reanalysed by Lee & Park (1999), with the inclusion of X-ray catalog 
data) . 

Owing also to the different definition of clusters in different observational 
samples or simulations, it is usual to work out from the cumulative cluster lu- 
minosity functions n{> L) the average distance between nearby objects Dl = 
n~^/^(> L) and to seek the dependence ofrc on D^. For 30 /i~^Mpc^ -Dl^ 60/i~^Mpc 
{h is the Hubble constant in units of 100 km/s/Mpc), the Bahcall & West 
(1988) conjecture that 

Tc ^ OADl (1.1) 

approaches most data, although, clearly, it cannot be considered an obser- 
vational output. However, over greater scales, there is no agreement among 
results obtained from different samples and it is not clear how much this may 
depend on the different cluster definitions. For Abell clusters, however, eq. 
(1.1) approaches data up to ~ 90/i~^Mpc. In this work, we shall assume the 
validity of eq. (1.1), as an example of gross discrepancy from assuming that 
Tc is L-independent. 

In the presence of such kind of L-dependence, two questions arise: (i) As a 
byproduct of the Limber equation, a scaling law holds, among the coefficients 
of the angular functions for samples with different limiting magnitudes. Is 
such scaling law still valid? We shall show that the same scaling law holds 
in the case (1.1), as in the constant Tc case, (ii) If the Limber equation is 
then formally used to work out the spatial clustering length, which value 
of Tc is obtained? We shall show that a naive use of the standard Limber 
equation leads to working out an apparent correlation length, close to the 
average distance between nearby clusters for the whole sample. 

As is obvious, our treatment applies to any kind of objects, not to clusters only. 
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Furthermore, the law (1.1) can be relaxed in various ways, obtaining related 
results, provided that the luminosity dependence of the clustering length (r^,) 
can be expressed through a function of the mutual distance of the objects be- 
longing to an intrinsical luminosity limited set (-Dl). Such generalizations will 
not be debated here. Our treatment also neglects departures from euclidean 
geometry and K-correction. Their inclusion does not change the conclusions 
of the paper, although slightly complicating final expressions; of course, for 
several detailed applications, they are to be included, but, at this stage, they 
are unessential to make our points. 

In order to show our points, we shall first need to reobtain the Limber equation 
in a way slightly different from textbooks (see, e.g., Peebles 1980a); this is the 
content of section 2 and 3. Using algorithms and expressions deduced in these 
sections, in section 4 we show the points made in the Introduction. Section 5 
contains a few final remarks. 



2 Sample definition 



The Limber equation was introduced to work out the spatial 2-point corre- 
lation function from angular data assuming that: (i) the 2-point function is 
luminosity independent; (ii) there exists a universal luminosity function for 
the class of objects considered. 

If the objects are galaxies, the luminosity function can be given the Schechter 
form (Schechter 1976): 



Here 1.1-1.2) is a phaenomenological parameter; L^, is a typical galaxy 
luminosity, close to the top luminosity of the sample. The number density of 
galaxies with luminosity exceeding L reads 

oo 

n{> L)^ j dL(t){L) = (t),U{L/L,) (2.2) 

L 

with 

oo 

[/(A) = j duu-"" ei^^{u) (2.3) 

A 

and, henceforth, 0(L) = —dn{> L)/dL. 

If we deal with clusters or other kinds of objects, the expression (2.1) may 
no longer be a fair approximation. In the sequel, however, we shall only need 
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that n{> L) can be expressed as in eq. (2.2), even though U{\) has not the 
expression (2.3). In eq. (2.2), (f)^ should be expressed in Mpcr^h^ and approach 
the total number density for the class of objects considered and, of course, 
there must also exist a typical luminosity L*, even though the distribution 
around it is not given by eq. (2.1). As previously stated, let us then define 

Dl = n-'/%> L) ■ n{> L) = D^^ , (2.4) 

clearly dn/dD = 

The apparent luminosity of a source of intrinsical luminosity L, at a distance 
r from the observer, is / = L/inr^ and the depth d^ of a sample of objects, 
with apparent luminosity />/,„, is defined so that 

Im = L^/^ndl . (2.5) 

Hence, requiring that an object at distance r has / > Im, means that its 
intrinsical luminosity 

L > U{r/d,f . (2.6) 
For instance, the angular density of objects, in a sample of depth d*, reads 

<x> oo oo 

nn(d.) = / drr' J dL<t>{L) = j drr'Dl^^^^,^^, = dlG , (2.7) 

L»{r/d*y 

where 

oo 1 °° 

G = I dqq'Dll^, = ^TsT^ / dL^Dl^ (2.8) 



is a universal constant, provided that (p{L) is a universal distribution, quite 
independently from the exphcit form of U{\). According to eq. (2.7), uq (x dl, 
a well known volume effect. 



3 The case of luminosity independent correlations 

All above relations hold provided that the assumption (ii) is true. For a sample 
of objects of depth d^, the angular 2-point function 

oo oo oo oo 

w(e) = n^^ I dnrl I dL,(f>(L,) J dr2rl j dL20(L2) C(ri2, L^, L^) 

L.(r-i/d*)2 U{rild»f 

(3.1) 

is obtainable from the spatial 2-point functions ^(ri2, i^i, ^^2) and from the 
luminosity function 0(L). 
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Let us first review the standard case, when also the assumption (i) holds, 
i.e. C = (ro/r)'^ with constant (luminosity independent) 7 and Tq (correlation 
length). In this case, changing integration variables, eq. (3.1) immediately 
yields: 



00 00 ^00 



w{e) = n^\l j dnrl J dD.^ Jdr^rl J dD.^n,^ (3.2) 



and, therefore, 

HO) = n^'r2 1 / ^^2r,^i?4,^/,,).r,7 (3.3) 



Herefrom, using eq. (2.7) and performing the changes of variables: 

2r = ri + r2 , q^r/d^ , ud^O = ri - r2 (3.4) 
in order that ri2 — ud^9{q^ + u'^y/^, we work out 

^^) = A(^)'^'-" (3.5) 



with 
where 



+00 



y du{l + uY^/^ (3.7) 
—00 

is a purely numerical constant. 

In above relations, the distances Dl were used, instead of the number densities 
n{> L) — D^^, which might even seem more convenient, in order to facilitate 
the forthcoming passage to the case of luminosity dependent correlations. 

Eq. (3.5) puts in evidence the scaling properties of the angular 2-point func- 
tion. The only dependence on the depth of the sample (d*) is factorized, as 
is a universal expression (still, provided that (p{L) is a universal distribu- 
tion). It is often stated that such scaling is arises from the fact that Tq does 
not depend on luminosities. Here we shall verify that such scaling is a more 
general property. 
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4 The case of luminosity-dependent correlations 



In the case of galaxy clusters, we shall now assume that the Bahcall and West 
conjecture holds: 

e(ri2,>L)=a(L>i/ri2r . (4.1) 

Here the numerical constant a ~ 0.4 and D^^ — n{> L) is the number density 
of clusters with luminosity > L (see eq. 2.5). Eq. (4.1) is an integral corre- 
lation law, as it concerns all objects with luminosities above L. In principle, 
it is possible to define also a differential correlation law, holding if objects of 
intrinsical luminosities Li and L2 only are considered. Such differential law is 
not to be measured, but we want to show that, the integral expression (4.1) 
follows from assuming that the differential law 

C(ri2, Li, L2) = U4D^2/ri2V (4.2) 

holds, provided that a = C*(l — 7/6)^^. In fact 

00 00 
e(ri2, >L) = n-\> L) ^ J dL^<P{L,)Dl{' J dL2ct>{L2)D]!^\l^ (4.3) 



and, changing variables as we did to pass from eqs. (3.1) to (3.2), 

\ 
I 



dD 



2)4-7/2 



2 

-2 / 



C.(l-7/6)-M — ) (4.4) 

T 12 



according to eq. (4.1). 



Let us now use the expression (4.2) to work out the angular function. In this 
case the correlation length has a precise dependence on luminosities. In spite 
of that, we shall find that eq. (3.5) is still formally true and, in particular w{9) 
has the same scaling properties as in the case of luminosity independent tq. 

In fact, if we replace the expression (4.2) in eq. (3.1) and perform the opera- 
tions leading there to eqs. (3.2) and (3.3), we find 

00 00 3 °° 3 

w{e) = n^^C. j drirl j dD^^^ j dr2rl j dD^^^r^^ 



00 00 

nn'W - 7/6)-^ / c^rir?D-JX. / dr^rlD-^^.rg (4.5) 
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Performing here again the changes of variables (3.4), we obtain 

oo +00 

w(e) = n^' C(l - l/6)-'d're'-^ I dqq'Dl'J;^ J du (q' + u^)'^!^ (4.6) 

-00 

and, using eq. (2.7), we return to 



provided that 



I r ^rvx - ^ Q ^ ^ ^ Cry -/-/j, ^ , o 



6+7 



^74 — r " ^"^ - 9r 7- 7/2 JO (Ar.^ 



(c-y is defined in eq. 3.7). 

Besides of finding the same scahng properties as in the case of L-independent 
To, we can now try to answer a practical question: If we interpret angular 
results, possibly because of the scaling w oc dj, as originating from a L- 
independent Tq and work out the Tq value, what value is obtained? 

This question is answered by requiring that 

r^,A, = f^A, (4.9) 

and, using eqs. (3.6) and (4.8), this yields 



^0^- .oo...._w.^-6 . (4-10) 



which is the apparent value of tq. 

Eq. (4.10) can be used either with simulations, working out from them and 
performing numerically the two integrations, or with analytical expressions. In 
the case of the Schechter law (2.1), using the function U (A) defined in eq. (2.3), 
it is easy to see that 



/0°°rfAA2-^/2[/2-7/3(A) 



1/7 



/o~dAA2-7/2[/2(A) 
The expression in square brackets can be turned into: 

,n/3 , , dx exp(-x) X- {I- dt eM-t) ^7(2^)'"'^' 



(4.11) 



/o°° dx exp(-x) X-- (/o~ dt eM-t) 
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7 


Cl 


C2 


1.6 


-1.95 


3.29 


1.8 


-1.93 


3.25 


2.0 


-1.91 


3.20 


2.2 


-1.88 


3.15 



with ei = 2 - 7/2 + (2 - 7/3)(l - a) and 63 = 2 - 7/2 + 2(1 - a), while 
v(q;,7) = 2'^2"'"-^/(2 — 7/3)'^^"'"^. Using a Gauss-Laguerre integration algorithm, 
double integrals can be then turned into double summations and analytically 
evaluated. For the narrow a interval, for which such expression is needed, 
however, we can verify that the linear expression 

ro/^y^ = Cl *q; + C2 (4.12) 

is better approximated than 0.2 %. 

In Tabic 1 we give ci and C2 for the relevant 7 interval. Using such values 
it is easy to see that the discrepancy between the correlation length ro and 
the average distance between nearby objects is expected to be ^ 10%. Such 
result, clearly, depends on the use of the Shechter law (2.1), which also allows 
a confortable numerical integration. We have considered a few other expres- 
sions, e.g. the expression deduced from the Press & Schecter (PS) cluster mass 
function, assuming a constant spectral index rit for the transfered fluctuation 
spectrum and L oc M^. In most such cases, numerical integration only is 
possible and integration details are cumbersome. For the PS case, however, 
the expression (4.12) give results approximated up to some percents, provided 
that we assume 0^ = 27r~^/^nti/*, where u* = p^/M^ (p^n is the average matter 
density and is the mass attributed to the cluster with typical luminosity 



5 Final remarks 

According to Peebles (1980a), the scaling relation (3.5) played an important 
role in testing that the angular correlations of galaxies in the catalogs do reflect 
the presence of a uniform spatial galaxy clustering, rather than something else, 
e.g. systematic errors due to patchy obscuration in the Milky Way. This work 
does not contradict such statement, corresponding to the point (ii) at the 
beginning of section 2. However, Hauser and Peebles (1973), in their seminal 
work on cluster clustering, made use of the scaling relation (3.5), as a test for 
the cluster correlation length Tc ~ 30 h~^Mpc they obtained for Abell clusters. 
This test, as we showed, has only a partial validity. 
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For the sake of curiosity, let us outline that, if the galaxy luminosity function 
values, 0* ~ 10~^ and a ~ 1.2 are used in eq. (4.11), we obtain a galaxy 
correlation length ~ 6 /i~^Mpc, for 7 values in the interal 1.6-1.9. This 
output should not be overstressed, as it has been known for a long time that 
the average distance between nearby galaxies Dg and the galaxy correlation 
length Tg have close values. In the case of galaxies, there is little evidence of 
a linear relation between values and Dg values for intrinsical magnitude 
limited samples. However, as galaxy luminosities essentially arise from stellar 
populations, they may be subject to significant variations, as their stellar 
and dust contents evolve on time scales certainly longer than observational 
times (see, e.g., Silva et al 1999). It is possible that the proximity of Vg and 
Dg indicates that, averaging luminosities over suitable time scales, a relation 
similar to clusters holds also for galaxies. 

In the case of galaxy clusters, Dalton et al (1997) investigated the scaling 
properties of the angular functions for two different samples with apparent 
magnitudes m > 18.94 and 19.3 < m < 19.5, for which they had obtained 
suitable U{X) expressions. Assuming a constant spatial clustering length Tc, 
they obtained that the predicted scaling properties hold, while 7 = 2.14 and 
Tc = 14.3 h~^Mpc. According to this work, finding a fair scaling does not 
require a constant Vc and, henceforth, any conclusion based on its existence is 
somehow premature. 
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